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FINAL  TECHNICAL  REPORT 


This  is  the  final  report  for  AFOSR  Grant  Number  F49620-93-1-0186,  a  three-year  grant 
to  the  University  of  Maryland  that  started  on  March  1,  1993  and  concluded,  after  a  no-cost 
extension,  on  June  15,  1996.  The  grant,  entitled  “Active  Control  of  Surge  and  Stall  in  Axial 
Flow  Compressors,”  focused  on  modeling  and  control  design  issues  for  the  nonlinear  surge  and 
rotating  stall  phenomena  that  can  occur  in  axial  flow  compressors  operated  in  the  regimes  of 
highest  performance.  The  cognizant  AFOSR  Program  Manager  for  the  grant  was  Dr.  Marc 
Q.  Jacobs,  of  the  Directorate  of  Mathematical  and  Information  Sciences  and  Computational 
Mathematics.  The  Principal  Investigator  on  the  grant  was  Dr.  Eyad  H.  Abed,  a  Professor  in 
the  Department  of  Electrical  Engineering  who  holds  a  joint  appointment  with  the  Institute 
for  Systems  Research.  Dr.  Raymond  A.  Adomaitis,  an  Assistant  Research  Scientist  with  the 
Institute  for  Systems  Research,  was  Co-  Principal  Investigator.  (Recently  Dr.  Adomaitis 
accepted  an  appointment  as  Assistant  Professor  of  Chemical  Engineering  at  the  University 
of  Maryland.) 

The  technical  report  begins  with  a  summary  of  the  main  goals  of  the  project.  Next,  the 
main  achievements  are  discussed.  The  remainder  of  the  technical  report  gives  information 
including  personnel  supported  under  the  grant,  publications  benefiting  from  the  grant,  in¬ 
teractions  and  transitions,  and  honors  and  awards  received  by  the  co-PIs  during  the  term  of 
the  grant. 

1  Overview  of  the  Research  Effort 

The  goals  of  the  research  effort  may  be  broadly  divided  into  five  categories,  as  follows.  Note 
that  the  first  four  of  these  categories  were  in  the  original  proposal,  while  the  fifth  was  added 
during  the  second  year  of  the  effort.  This  fifth  category  was  found  to  be  beneficial  to  the 
control  design  aspects  of  the  work. 

(i)  Nonlinear  Dynamics  of  Compression  Systems.  This  aspect  of  the  research  includes 
analytical  and  numerical  study  of  the  qualitative  changes  in  the  dynamical  behavior 
of  compression  systems  that  occur  as  parameters  are  varied.  This  results  in  an  under¬ 
standing  of  both  the  instability  mechanisms  responsible  for  stall  and  surge  inception, 
and  the  nonlinear  growth  to  fully-developed  surge  and  rotating  stall. 

(ii)  Synthesis  of  Robust  Nonlinear  Controllers.  This  aspect  of  the  research  mainly  involves 
application  of  bifurcation  control  techniques  to  active  stall  control  design,  allowing 
operation  near  peak  pressure  rise.  It  also  involves  extensions  of  bifurcation  control 
results  as  needed. 

(iii)  Global  Stability  Analysis  of  Controllers.  This  goal  of  the  research  effort  involves  non¬ 
linear  stability  analysis  of  models  of  controlled  compression  systems. 

(iv)  Modeling  Enhancements.  This  phase  of  the  research  involves  development  of  high- 
fidelity,  multi-stage  compression  system  models  along  with  scalable  solution  algorithms 
suitable  for  massively  parallel  processors.  Emphasis  is  placed  on  modeling  techniques 
which  provide  a  rational  framework  for  generating  reduced-order  models  useful  for 
bifurcation  analysis  and  control. 


(v)  Modal  Participation  and  Instability  Precursors.  This  phase  of  the  research  involves 
two  related  aspects:  study  of  the  participation  of  states  in  modes  of  instability,  and 
algorithms  for  real-time  detection  of  impending  instability.  In  addition,  it  is  impor¬ 
tant  to  apply  these  results  to  multistage  compression  system  models,  and  to  deduce 
implications  for  sensor  and  actuator  placement  in  multistage  compression  systems. 

2  Accomplishments/New  Findings 

The  main  accomplishments  of  the  research  project  can  be  divided  into  results  on  modeling 
and  and  results  on  control  of  multistage  compression  systems.  Moreover,  the  modeling 
results  and  the  control  design  results  both  can  be  further  viewed  as  numerical  or  analytical. 
Thus,  detailed  numerical  models  as  well  as  less  detailed  analytical  models  were  developed. 
The  numerical  models  developed  also  include  graphical  display  capabilities  for  illustration 
of  the  evolution  of  the  gas  flow  fleld  in  a  compression  system.  There  are  important  and 
immediate  benefits  to  the  parallel  development  of  numerical  and  analytical  approaches  for 
the  analysis  and  control  of  compression  systems.  For  example,  it  is  very  difficult  to  design  a 
controller  based  on  a  numerical  model  of  the  dynamical  system  to  be  controlled,  especially  if 
the  system  is  complex.  However,  availability  of  a  less  accurate  analytical  model  that  reflects 
the  structure  of  the  system  and  any  fundamental  nonlinearities  can  facilitate  the  design  of  a 
controller.  The  controller  designed  in  this  way  can  then  be  tested  on  the  numerical  model. 
In  fact,  the  numerical  model  can  be  used  to  “tune”  the  control  gains.  This  methodology  can 
also  be  applied  to  an  experimental  set-up  rather  than  a  numerical  model. 

In  the  following,  we  list  the  main  accomplishments  of  the  research  project. 

1)  We  have  developed  a  modeling  approach  unique  in  the  field  of  turbocompressor  simulation 
in  that  not  only  do  the  forcing  functions  depend  on  spatial  position,  but  they  are  also 
computed  solely  from  the  equations  of  fluid  motion  and  the  geometry  of  the  compressor 
blading.  The  simulation  model  developed  was  employed  to  study  hysteresis  in  the  compressor 
behavior.  An  important  aspect  of  this  work  is  studying  the  predictive  capabilities  of  such 
a  model  in  simulating  the  inception  of  rotating  stall  as  the  stall  margin  is  crossed,  and 
following  the  growth  of  these  instabilities  to  fully-developed  rotating  stall.  This  work  has 
shown  that  simulations  of  rotor  blade  wake  momentum  defect  growth  to  a  fully  separated 
flow  field  blocking  a  blade  passage  is  consistent  with  the  conceptual  model  of  compressor 
stall,  originally  proposed  by  lura  and  Rannie  (Trans.  ASME,  1954).  Stall  cell  propagation 
rates  of  approximately  1/2  rotor  speed  were  observed  in  our  simulations.  The  stall  cell 
propagation  mechanism  put  forth  by  the  cited  study  is  corroborated  by  our  simulations.  An 
important  step  towards  understanding  turbocompressor  flow  instability  inception  and  control 
that  has  been  completed  in  this  work  is  the  development  of  an  accurate  fluid  dynamic  model. 
Issues  such  as  identification  of  the  primary  instability  mechanism,  determining  whether  the 
onset  of  instabilities  can  be  detected,  and  comparisons  of  different  control  approaches  (one¬ 
dimensional  vs.  “distributed”  controls)  can  be  resolved  with  models  of  sufficient  fidelity. 

Because  of  the  complexity  of  multi-stage  compression  systems,  direct  modeling  of  flows 
at  the  individual  rotor/stator  blade-scale  level  is  computationally  impractical.  Furthermore, 
dynamic  models  of  such  detail,  even  if  computationally  tractable,  would  not  be  amenable 
to  theoretical  analysis  techniques  used  to  gain  insight  into  stall  inception  mechanisms,  and 
certainly  are  not  suitable  for  controller  development. 

Because  of  the  limited  predictive  capabilities  of  compressor  flow  simulators  in  the  past. 


a  large  body  of  experimental  data  on  compression  system  performance  has  accumulated. 
Theoretical  understanding  of  the  bifurcations  responsible  for  compressor  stall  and  surge  has 
evolved  in  parallel  with  this  experimental  work.  Our  modeling  work  is  meant  to  bridge  the 
gap  between  the  two,  incorporating  experimentally-determined  correlations  describing  flow 
features  too  difficult  to  simulate  from  first  principles,  combined  with  rigorous  flow  modeling 
where  possible.  In  the  numerical  modeling  work,  a  consistent  path  to  reduced-order  models 
was  always  kept  clear. 

Computing  a  representative  Reynolds  number  based  the  hub  diameter  length  scale  in¬ 
dicates  that  the  flow  through  the  compressor  is  highly  turbulent  This  means  compressor 
instabilities  are  transitions  in  the  structure  of  the  time-smoothed  flow  field.  With  this  in 
mind,  we  see  that  a  rational  model  of  compressor  flows  can  be  based  on  averaging  the  Navier- 
Stokes  equations  over  a  length  scale  longer  than  the  turbulent  fluctuations  and  shorter  than 
the  large-scale  structures.  The  dynamics  of  large-scale  structures  depend  on  compression 
system  geometry  and  other  specific  characteristics;  the  small-scale,  energy  dissipating  struc¬ 
tures  are  more  universal  and  so  closure  of  the  averaged  Navier-Stokes  is  provided  by  a  subgrid 
modeling  technique. 

To  overcome  the  computational  infeasibility  of  detailed  simulation  of  the  fluid  flow  dy¬ 
namics  at  individual  blade  levels  in  multistage  compressors,  we  approximate  the  power  input 
to  the  fluid  flow  field  in  each  stage  by  forcing  functions  to  the  fluid  equations  of  motion. 
These  data  are  available  in  the  form  of  experimentally-determined  stage-by-stage  perfor¬ 
mance  curves.  This  modeling  approach  forms  the  basis  of  the  joint  experimental/simulation 
research  outlined  in  the  Cooperative  Research  and  Development  Agreement  between  the  Uni¬ 
versity  of  Maryland  and  Army  Research  Office  Vehicle  Propulsion  Directorate.  The  goals  of 
this  project  are  to  understand  flow  instabilities  in  the  seven  axial,  single  centrifugal  stage 
T-55  compressor  and  development  of  a  stability  augmentation  system. 

In  the  CFD-based  modeling  discussed  above,  the  power  input  forcing  functions  were 
defined  for  each  rotor  and  stator  row.  If  we  spread  the  forcing  functions  over  the  entire  com¬ 
pressor  section  and  make  further  simplifying  assumptions  on  the  turbulence  modeling  and 
size  of  fiow  perturbations,  we  can  reduce  the  model  to  a  single  partial  differential  equation. 
Analysis  of  this  reduced  system  by  linearized  stability  analysis  and  weakly  nonlinear  meth¬ 
ods  is  possible.  The  findings  of  this  analysis  demonstrate  that  a  sequence  of  bifurcations 
corresponding  to  stalled  flows  of  different  cell  numbers  takes  place  along  the  uniform  flow 
solution  branch.  The  first  of  these  stall  inception  instabilities  is  a  sub  critical  bifurcation. 
Secondary  bifurcations  render  this  flow  pattern  stable,  leading  to  a  hysteresis  loop. 

2)  The  simulation  model  developed  was  employed  to  study  hysteresis  in  the  compressor 
behavior.  In  the  model,  the  throttle  was  closed  through  the  range  where  the  flow  was 
found  to  become  unstable.  Simulation  showed  the  onset  of  stall  soon  afterwards.  The 
disturbance  was  found  to  grow  to  a  large-amplitude  rotating  stall  cell  after  only  a  few  rotor 
revolutions,  growing  only  out  of  the  natural  flow  asymmetry  induced  by  the  rotor-stator 
blade  interactions.  After  a  period  of  disturbance  growth  to  fully  developed  rotating  stall, 
the  throttle  opening  was  returned  to  its  original  size,  which  had  initially  corresponded  to  a 
locally  asymptotically  stable  uniform  flow  operating  point.  After  reaching  equilibrium,  we 
found  that  the  rotating  stall  persisted,  a  clear  indicator  of  multistability  (hysteresis). 

3)  At  the  inception  of  instability  or  stall  in  a  multistage  compression  system,  one  would 
expect  that  not  all  stages  participate  equally  in  the  instability.  It  is  important  to  determine 
the  stage  or  stages  that  are  most  impacted  by  the  instability.  One  main  reason  for  this 


is  that  such  knowledge  is  useful  for  actuator  placement.  Actuators  are  best  placed  in  the 
physical  neighborhood  of  an  instability  onset. 

We  have  applied  Selective  Modal  Analysis  (SMA)  to  this  issue.  SMA  was  originally 
introduced  in  the  electric  power  systems  field,  where  it  has  been  useful  in  determination  of 
the  physical  nature  of  modes  of  instability  (e.g.,  electrical  vs.  mechanical).  It  has  also  proved 
useful  in  order  reduction,  and  in  sensor  and  actuator  placement  issues.  In  this  technique, 
each  mode  of  a  dynamic  system  is  related  to  subsystem  modes,  and  thus  to  state  variables 
of  the  system.  The  method  results  in  numerical  indices,  called  participation  factors,  that 
indicate  the  degree  to  which  any  state  variable  participates  in  a  system  mode.  This  feature 
allows  determination  of  the  participation  of  each  stage  in  an  overall  compression  system 
mode.  The  main  application  is  to  determination  of  the  stage(s)  most  involved  in  a  mode 
of  instability,  with  the  associated  implications  for  control  design.  Particularly  important 
among  these  are  implications  for  actuator  placement.  For  example,  insights  can  be  obtained 
into  which  bleeds  or  variable  stator  blade  rows  to  actuate. 

4)  Besides  applying  SMA  (in  particular  participation  factors)  to  multistage  compression 
system  models,  we  have  also  re-examined  the  concept  of  participation  factors  and  given  a 
derivation  of  the  expression  for  these  quantities  based  on  a  new  first-principles  definition. 
This  is  useful  not  only  from  a  conceptual  perspective,  but  also  as  a  basis  for  extending  this 
useful  method  to  nonlinear  problems  and  in  other  useful  directions. 

5)  We  have  introduced  the  use  of  “precursor  signals”  into  the  study  of  multistage  com¬ 
pression  system  instability.  Precursors  were  originally  studied  by  Weisenfeld  as  a  tool  for 
understanding  the  amplification  of  noise  in  systems  near  bifurcation  points.  We  have  ex¬ 
panded  this  notion  to  address  system  monitoring  and  signal  processing  techniques  that  result 
in  a  “warning  signal”  that  appears  when  system  parameters  drift  dangerously  close  to  critical 
(bifurcation)  values.  Past  these  critical  values  of  the  system  parameters,  the  system  oper¬ 
ating  condition  will  likely  lose  stability,  corresponding  to  stall.  The  availability  of  a  reliable 
precursor  signal  is  thus  very  important  in  the  safe  operation  of  a  physical  system  the  model 
for  which  is  highly  uncertain.  Further  work  on  precursors  usin  higher  order  spectral  analysis 
also  was  undertaken  and  this  work  continues  in  a  new  grant  from  AFOSR. 

6)  We  have  coupled  the  work  on  precursors  with  the  SMA  approach  mentioned  above,  result¬ 
ing  in  a  method  for  deciding  the  best  location  for  placing  sensors  for  continuous  monitoring 
of  a  multistage  compression  system.  Briefly,  it  was  found  that  compression  system  stages 
whose  variables  have  high  participation  factors  provide  significantly  larger  precursor  signals 
as  compared  to  other  stages. 

7)  We  have  continued  our  investigations  of  control  of  systems  with  bifurcations  with  appli¬ 
cation  to  axial  flow  compressors.  Having  shown  that  stall  inception  takes  place  by  means  of 
a  subcritical  bifurcation  near  peak  pressure  rise,  we  can  define  the  control  problem  as  one  of 
reducing  the  hysteresis  loop  to  improve  operability.  This  type  of  control  effectively  reduces 
the  “nonlinear”  stall  margin  defined  by  the  range  of  multistability.  Using  the  analytical  mul¬ 
tistage  models  and  the  numerical  reduced-order  models,  we  have  developed  stall  suppression 
methods  by  direct  control  of  the  instability  inception  bifurcations.  Further  refinements  of 
the  stall  controllers,  such  as  the  benefits  of  nonsmooth  controls,  have  been  investigated. 

We  have  also  analyzed  controller-induced  oscillations  found  in  experiments  at  United 
Technologies  Research  Center  (UTRC).  These  oscillations  occur  when  a  basic  bifurcation 
control  law  proposed  by  our  group  is  applied  in  high  speed  compressors.  The  UTRC  re¬ 
searchers  then  modified  the  control  law  to  alleviate  these  oscillations.  We  have  studied  the 


reasons  for  the  occurrence  of  the  oscillatory  behavior  and  means  for  alleviating  it  analyti¬ 
cally.  This  analytical  study  makes  use  of  tools  from  bifurcation  theory  and  singular  pertur¬ 
bation  theory.  Briefly,  we  have  shown  that  as  the  compressor  speed  parameter  (Greitzer’s 
B-parameter)  is  increased,  the  system  differential  equations  become  singularly  perturbed 
and  a  Hopf  bifurcation  that  occurs  on  the  bifurcated  equilibrium  solution  moves  arbitrarily 
close  to  the  original  bifurcation  point. 

We  have  also  developed  a  control  design  approach  using  the  harmonic  balance  method. 
We  have  given  explicit  calculations  for  approximate  determination  of  Hopf  and  period  dou¬ 
bling  bifurcation  points  in  nonlinear  systems,  and  for  determination  of  the  effect  of  control 
action  on  these  bifurcations.  Other  papers  in  the  list  of  publications  address  further  appli¬ 
cations  of  bifurcation  control  to  convection  and  power  systems  and  to  the  control  of  chaos. 

8)  We  have  contributed  to  the  limited  and  relatively  new  literature  on  the  control  of  systems 
with  bifurcations  by  writing  several  review  articles  on  the  subject,  one  of  which  appears  in 
the  recently  published  Control  Systems  Handbook  from  CRC  Press.  These  articles,  which 
are  written  for  a  general  control  systems  audience,  are  expected  to  help  increase  the  under¬ 
standing  of  the  role  of  bifurcations  and  nonlinear  dynamics  in  control  system  analysis  and 
design.  We  have  also  prepared  a  set  of  lecture  notes  on  the  subject  that  has  been  used  in  a 
graduate  special  topics  course  at  the  University  of  Maryland  as  well  as  at  special  workshops 
in  Italy  (University  of  Florence)  and  Taiwan  (National  Chiao-Tung  University).  These  notes 
are  being  revised  to  prepare  them  for  publication  in  book  form. 
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5  Interactions/Transitions 

1)  In  the  final  year  of  the  grant,  an  AFOSR  “Monthly  Miracle”  article  was  prepared  in 
cooperation  with  the  Program  Manager  (Dr.  Marc  Jacobs)  and  with  Dr.  Carl  Nett  of  UTRC. 
The  article  later  appeared  in  the  annual  publication  “AFOSR  Research  Highlights.”  The 
article  describes  the  transition  of  a  compressor  control  concept  from  this  group  to  industry. 

2)  An  invited  talk  entitled  “Intelligent  Control  Issues  for  PEBB-Based  Systems”  was  given 
by  E.H.  Abed  at  the  ONR-Virginia  Tech  Control  Technology  Workshop:  The  Regulation 
and  Distribution  of  Power  in  Large  Systems,  held  at  the  Virginia  Center  for  Innovative 
Technology,  Herndon,  VA,  24-25  April  1995.  The  nonlinear  dynamics  issues  that  may  arise 
in  networks  of  PEBBs  (power  electronic  building  blocks)  that  ONR  has  set  as  a  research 
goal  were  discussed.  Some  of  the  lessons  learned  in  the  current  grant  were  communicated. 


3)  A  talk  entitled  “Precursors  for  nonlinear  instabilities,”  was  given  at  the  1995  AFOSR 
Grantees/ Contractors  Meeting  in  Dynamics  and  Control,  held  in  Minneapolis  during  5-7 
June  1995. 

4)  A  Cooperative  Research  and  Development  Agreement  with  the  Army  Vehicle  Propulsion 
Directorate  at  NASA’s  Lewis  Research  Center  is  in  place.  This  agreement  is  for  the  joint 
development  of  a  CFD  model  of  the  eight-stage  axi-centrifugal  compressor  for  the  T55  com¬ 
pressor,  found  in  the  turboshaft  engine  which  powers  the  Army’s  Chinook  helicopter.  This 
dynamic  model  will  be  used  in  developing  and  simulating  stall  controllers  for  the  T55. 

5)  There  has  been  significant  interaction  with  Dr.  Wishaa  Hosny  of  General  Electric  Aircraft 
Engines  on  multistage  compression  system  modeling.  Discussions  have  also  taken  place  with 
researchers  from  Wright-Patterson  Air  Force  Base  and  United  Technologies  Research  Cen¬ 
ter.  There  was  also  regular  communication  with  the  engineering  staff  of  Textron-Lycoming. 
Possible  uses  of  the  ideas  developed  in  this  project  for  control  of  combustion  instabilities  and 
for  turbine  modeling  and  control  are  being  discussed  with  industrial  researchers. 

6)  Some  of  the  ideas,  analysis  techniques,  and  simulation  packages  developed  in  the  context 
of  stall  control  research  carry  over  to  materials  manufacturing  process  control.  This  is 
a  growing  area  of  research  at  the  Institute  for  Systems  Research,  and  collaborative  work 
has  been  initiated  with  the  Materials  Directorate  at  WPAFB,  NRL,  North  Carolina  State 
University,  and  other  academic  and  government  laboratories. 

7)  Compressor  control  problems  have  been  discussed  with  faculty  of  other  universities  work¬ 
ing  in  this  area. 

8)  Dr.  Abed  co-organized  a  Workshop  on  Nonlinear  Control  and  Control  of  Chaos  that  was 
held  at  the  International  Centre  for  Theoretical  Physics  in  Trieste,  Italy,  during  June  17-28, 
1996.  The  workshop  brought  together  leading  researchers  in  the  area  of  control  of  nonlinear 
dynamics  to  promote  cross-fertilization  of  ideas  from  the  engineering  and  science  disciplines. 
At  the  workshop.  Dr.  Abed  gave  a  lecture  entitled  “Control  of  Nonlinear  Dynamics  in 
Engineering  Systems.” 

6  Honors  /  Awards 

During  the  term  of  the  grant,  the  Principal  Investigator  (Eyad  H.  Abed)  was  the  recipient 
of  the  following  honors  and  awards: 

1.  The  O.  Hugo  Schuck  Best  Paper  Award  was  received  in  1994  by  Dr.  Eyad  H.  Abed. 
The  awarding  organization  was  the  American  Automatic  Control  Council. 

2.  Listed  in  the  biographical  collections  American  Men  and  Women  of  Science,  Who’s 
Who  in  Technology,  and  Who’s  Who  in  Science  and  Engineering. 


